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We consider preheating in the theory \\<f) A + \g 2 $ 2 X 2 1 where the classical oscillating inflaton 
field 4>(t) decays into x-P ar ticles and (^-particles. The parametric resonance which leads to particle 
production in this conformally invariant theory is described by the Lame equation. It significantly 
differs from the resonance in the theory with a quadratic potential. The structure of the resonance 
depends in a rather nontrivial way on the parameter g 2 /X. We construct the stability/instability 
chart in this theory for arbitrary g 2 /A. We give simple analytic solutions describing the resonance in 
the limiting cases g 2 /A <C 1 and g 2 /\ ^> 1, and in the theory with g 2 — 3A, and with g 2 = A. From 
the point of view of parametric resonance for the theories with g 2 = 3A and with g 2 — A have 
the same structure, respectively, as the theory |A</) 4 , and the theory j^(4> 2 ) 2 of an TV-component 
scalar field cj>i in the limit TV — *■ 00. We show that in some of the conformally invariant theories 
such as the simplest model |A(/> 4 , the resonance can be terminated by the backreaction of produced 
particles long before {x 2 } or (4> 2 ) become of the order <f> 2 . We analyze the changes in the theory of 
reheating in this model which appear if the inflaton field has a small mass. 
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I. INTRODUCTION 

The theory of reheating of the universe is one of the 
most important and least developed parts of inflationary 
cosmology. Recently it was found that in many realis- 
tic versions of chaotic inflation reheating begins with a 
stage of parametric resonance |J. At this stage the en- 
ergy is rapidly transferred from the inflaton field to other 
scalar and vector fields interacting with it. This process 
occurs far away from thermal equilibrium, and therefore 
we called it preheating. The theory of preheating is rather 
complicated. In [Q we gave only a brief summary of its 
basic features. A detailed investigation of preheating in 
the simplest chaotic inflation model describing a massive 
inflaton field <f> interacting with a massless scalar field \ 
was contained in our recent paper [^. It was found, in 
particular, that the resonance in such theories can be ef- 
ficient only if it is extremely broad. In such a situation 
preheating in an expanding universe looks like a stochas- 
tic process. 

In this paper we will concentrate on the theory of pre- 
heating in a class of conformally invariant theories such 
as j<p 4 + ^(px 2 ■ Different aspects of preheating in such 
theories have been studied in Refs. [l J^-|lO|. A specific 
feature of these models is that by a conformal transfor- 
mation one can reduce the investigation of preheating in 



these theories in an expanding universe to a much sim- 
pler theory of preheating in Minkowski space-time |l| . As 
a result, the parametric resonance does not exhibit the 
stochasticity found in M. However, stochastic resonance 
may appear again at the late stages of preheating if the 
fields <fi and x have bare masses which break conformal 
invariance. 

We will investigate preheating in the theories of the 
type of jcf) 4 + \^x 2 f° r various relations between the 
coupling constants g 2 and A. During this investigation 
(see specifically Sec. V and XII), we will discuss how 
the results of the previous papers on this subject are 
related to the picture which emerges from the current 
study. We will show that the development of the reso- 
nance in the various conformally invariant theories can 
be very different, depending on the particular values of 
parameters and the structure of the theory. For exam- 
ple, the model |0 4 -I- ^<\> 2 x 2 with g 2 = A or g 2 = 3A has 
only one instability band, but the structure of the bands 
and the characteristic exponents fik are completely dif- 
ferent. It is enough to change the ratio g 2 /\ only slightly, 
and the number of the instability bands immediately be- 
comes infinitely large. For this reason, it is dangerous 
to extrapolate the results obtained for a theory with one 
choice of parameters to a theory with another choice of 
parameters. As we will see, not only is the structure of 
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the resonances different in different models, but the self- 
consistent dynamical evolution with an account taken of 
the backreaction of produced particles can also be quali- 
tatively different. 

The main purpose of the present paper is to study the 
structure of the parametric resonance in the conformally 
invariant theories. These theories may describe many 
bose fields \ interacting with the inflaton field <fi with 
different coupling constants: 
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Here Xm stands for the m-th scalar field interacting with 
the inflaton field with the coupling constant g mi and in- 
teracting with curvature R with the coupling constant 
£ m . The equation for fluctuations in this general model 
( |i1f ) unifies the equations for fluctuations in the confor- 
mal models mentioned before including 5cj> fluctuations. 

Strictly speaking, this model is conformally invariant 
only for a specific choice of the parameters £ m : £ m = g. 
Nevertheless, in this paper we will consider the simplest 
models with £ m = 0. As we will see shortly, this dif- 
ference is not going to be very important because the 
average value of R vanishes when <f> <C M p . 

We will see that for the conformally invariant theories 
the only parameter actually responsible for the struc- 
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ture of the resonance for the field Xm is the ratio 
Furthermore, we will find that the strength of the reso- 
nance and the number and widths of the instability bands 
for the field \ m in the theory (0) depends on %»• non- 
monotonically. To get a general picture, we will construct 
the stability/instability chart for the equation for fluctu- 
ations on the two dimensional plane (fc 2 , ^-), see Fig. 3. 
The stability/instability chart gives us insight into the 
structure of the resonances in the conformally invariant 
theories. From this it will immediately be clear which of 
the fields Xn of (|l|) will be most amplified during preheat- 
ing. The stability/instability chart unifies our knowledge 
of the resonance for the various conformal models thus 
far considered in the literature. 

Note that the class of theories we are going to in- 
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vestigate include in particular the theory jC%2 4*1) of 

an iV-component scalar field </>,;. This theory has O(N) 
symmetry. One can identify the inflaton field <fi in this 
theory with the field <j>\. Then the quantum fluctua- 
tions of this field, just like the quantum fluctuations in 
the theory of a one-component field j4> 4 , will have effec- 
tive mass squared 3A0 2 , whereas the fluctuations of all 
other fields will have effective mass squared X(f> 2 . There- 
fore, the equation for the growth of the fluctuations of 
the field <fi = 4>\ (neglecting backreaction) will coincide 
with the equation for the growth of fluctuations of the 
field x coupled to the field </> with the coupling constant 



g 2 = 3A. Meanwhile, the equation for the growth of the 
fluctuations of the fields i ^ 1, will coincide with the 
equation for the growth of fluctuations of the field x with 
the coupling constant g 2 = A. This regime is especially 
important in the limit N — > oo, where the main contri- 
bution to particle production is given by the modes with 
i 1. Thus, the cases g 2 = A and g 2 = 3A are especially 
interesting and deserve careful investigation. 

This paper is organized as follows. In Sec. II we will 
describe the evolution of the background inflaton field 
4>(t) after inflation in the theory with the effective poten- 
tial V((j)) = \X(f> 4 . In Sec. II we will give its analytic so- 
lution. Then, in Sec. Ill, we derive the equations for fluc- 
tuations of the fields x and <f> in the conformally invari- 
ant theory, and reduce these to equations in Minkowski 
space-time. We show that these equations can ultimately 
be reduced to a single Lame equation with just one pa- 
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rameter, In Sec. IV we solve the resonance equations 

2 

numerically for an arbitrary and arbitrary momentum, 
fc, of fluctuations. This allows us to produce the main re- 
sult of our paper; we construct the stability/instability 
chart for fluctuations in the conformally invariant the- 
ories. In Sec. V we discuss the particular ranges and 
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values of the parameter ^- where the analytic methods 
for the description of the resonance can be developed. In 
Sec. VI - IX we perform an analytic investigation of the 
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resonance for some particular values of 4-. For different 

2 

values of ^- different analytic approaches will be devel- 
oped. We report a new method to treat the resonance 
when ^- = ; where n is an integer. We show that 

the solutions for 4- = "^" 2 +1 ' ) can be found in closed, ana- 
lytic form. This is done explicitly for the most interesting 
cases, n = 1 and n = 3 (i.e. for g 2 = A and g 2 = 3A), 
in Sec. VI, VII, and the Appendix. We also consider the 
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two opposite limits ^- <C 1 and ^- ^> 1 in Sec. VIII 
and IX respectively. Sec. X contains a discussion of the 
self-consistent dynamics of the system including backre- 
action of the created particles. In Sec. XI we describe 
the restructuring of the resonance which occurs when the 
backreaction is incorporated into the equations for fluc- 
tuations. We show that this is the leading effect which 
terminates the resonance in the theory j(f> 4 . In Sec. XII 
we discuss the modifications of the theory of preheating 
which appear when the inflaton field <j) is massive. This 
allows us to unify the results obtained in this paper with 
the results of our preceding investigation of preheating 
in the theory of the massive inflaton field Q . Finally, in 
Sec. XIII, we give a summary of our results and discuss 
their possible implications. 



II. EVOLUTION OF THE INFLATON FIELD 

We consider chaotic inflation with the potential 
V{4>) — iA(/) 4 . During inflation the leading contribution 
to the energy-momentum tensor is given by the inflaton 
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scalar field <f>. The evolution of the (flat) FRW universe 
is described by the Friedmann equation 
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In the theory with the potential V(4>) — jXcfi 4 , which 
we consider in this paper, the inflaton oscillations are 
not sinusoidal. The amplitude $ of the oscillations of 
the scalar field 4> m the limit t — > oo approaches the 
asymptotic regime 



where H = a/a. Let us note one more useful relationship 
between H(t) and <fi(t) which follows from the Einstein 
equations 



H = - 
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The equation for the classical field 4>(t) is 
4> + 3H <j) + Xct) 3 = . 



(3) 
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M p , the fric- 



For sufficiently large initial values of (p > 
tion term, 3H(b, in (Q) dominates over </> and the poten- 
tial term in (0) dominates over the kinetic term. This 
is the inflationary stage, where the universe expands 
quasiexponentially, a(t) — ao exp (/ dtH(t)). With a de- 
crease of the field <fi below M p , the "drag" term 2>H <j> 
gradually becomes less important and inflation termi- 
nates at <f> ~ Mp/2. After a short stage of fast rolling 
down, the inflaton field rapidly oscillates around the min- 
imum of V(4>) with the initial amplitude $o ~ 0.1M p . 
Although this value is below the magnitude needed for 
inflation, it is still very large. 

The character of the classical oscillations of the homo- 
geneous scalar field depends on the shape of its poten- 
tial V{4>). In Ref. |^] we considered the theory with the 
quadratic potential V((f>) = ^m<j) 2 . In that theory the 
fluctuations are harmonic, <f>(t) — $(t)sinmi, with the 



amplitude decreasing like 
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a- 3 / 2 . The 



scale factor at the stage of oscillations is a(t) aot 2 ' 3 , 
and the energy density of the inflaton field decreases in 
the same way as the energy density of nonrelativistic mat- 
ter cx a~ 3 . 




FIG. 1. Oscillations of the field <j> after inflation in the 
theory sSL . The value of the scalar field here and in all other 
figures in this paper is measured in units of M p , time is mea- 
sured in units of (V r \M p )~ 1 . 
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where TV is the number of oscillations after the end of 
inflation. 

To make calculations in this theory, and in particular, 
to find the form of the oscillations, it is convenient to 
make a conformal transformation of the space-time met- 
ric and the fields. For this we use the conformal time 



V 



and the conformal field, 



dt 
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tp = a<j>. (7) 

In the coordinates (77, x) the Klein-Gordon equation (^) 
for ip is 
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where ' stands for the derivative with respect to the con- 
formal time, The Friedmann equation (||) in these 
variables is 
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As one can see from Eq. (|S|) , the equation of motion for 
the field ip in the time variable 77 does not look exactly as 

the equation for the theory jtp 4 in Minkowski space. In 

j 2 

order to achieve it one would need to add the term j^R 
to the Lagrangian. However, this subtlety is not very 
important. First of all, soon after the end of inflation 
one has jc/) 4 3> j^R, and Xp 3 ^> —ip. Moreover, it is 
known that the energy- momentum tensor of the field <fi in 
the theory j<p 4 , when averaged over several oscillations, 
is traceless (p — p/3) |fj"l|| . In this case one has R = 0, 
a(rj) ~ 77, and a" = 0, so that the last term in Eq. (||) 
vanishes: 



Xip 3 = 



(10) 



The Friedmann equation (^) averaged over several os- 
cillations of the field (j> in the regime (j> <C M p also takes 
a very simple form: 
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where we have introduced the conformal energy density, 
Pv = \p n + \p 4 . 



3 



It is convenient to express in terms of the amplitude 
<p of the oscillations of the field ip: p v = \<p> A . Equation 
( |To| ) has an oscillatory solution with a constant amplitude 
and the conformal energy p v . Then from (Jll|) we find 



0(77) 



-2 



2n\ (p 
~3~M 



V 



p 



-2 



7rA ft 

6~M 



(12) 



-</? m 



the 



As we expected, in this regime the last term 
equation_([lO]) vanishes. 

Eq. (|l(j) can be reduced to the canonical equation for 
an elliptic function. Indeed, let us use a dimensionless 
conformal time variable 
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Then we can rescale the function ip = acj) = ipf(x). The 
function f(x) has an amplitude equal to unity and obeys 
the canonical equation for the elliptic function. The in- 
tegral of this equation, f' 2 = |(1 — has the solution 
in terms of an elliptic cosine 
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As claimed, oscillations in this theory are not sinusoidal 
but are given by an elliptic function. The energy density 
of the field 4> decreases in the same way as the density of 
radiation, i.e. as et~ 4 . 



over a period is |. The potential energy density jXip 4 
averaged over a period of oscillation is equal to \pip, and 
the average kinetic energy \<p' 2 is given by \p v - 
The elliptic cosine can be represented as follows: 
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The amplitude of the first term in this sum is 0.9550; the 
amplitude of the second term is much smaller, 0.04305. 
The full solution ( |l4| ) is plotted in Fig. 1 (solid curve), 
alongside the leading harmonic term in the series jl^), 
cos0.8472x (dotted curve). Although the first harmonic 
term is very close to the actual form of oscillations, it will 
be important for the investigation of the general structure 
of stability/instability bands in this theory that f(x) is 
not exactly equal to cos 



III. EQUATIONS FOR QUANTUM 
FLUCTUATIONS OF THE FIELDS <j) AND X 

We will consider here the interaction between the clas- 
sical inflaton field, 0, and the massless, quantum scalar 
field, Xi with the Lagrangian (§). The Heisenberg repre- 
sentation of the quantum scalar field \ is 
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FIG. 2. The exact solution ( p^ ) for the oscillations of 
the inflaton field after inflation in the conformally invariant 
theory j\<f) 4 . We show the field in rescaled conformal field 
and time variables, f(x) — % (solid curve) and the first term, 
cos(0. 8472a;), in its harmonic expansion ( p^5[ ) (dotted curve). 

The solution ( p^ ) has some interesting properties 
which are not usually elucidated in the literature. It 
matches the solution describing the slow rolling of the 
field <p at the end of inflation if one takes xo « 2.44. The 
period of the oscillations (in units of x) is T = 4K(-^) = 
r 2 (i/4) 



7.416, K stands for the complete elliptic inte- 
gral of the first kind. The effective frequency of oscil- 
lations is 2tt/T « 0.8472 f|. The value of / 4 averaged 



where a,k and are the annihilation and creation oper- 
ators. For a flat Friedmann background with scale factor 
a(t), the temporal part of the eigenfunction with comov- 
ing momentum k obeys the following equation: 



a / k 2 
Xk + 3-Xk+ — + g 2 d? 2 )x,k = 
a \ a z 



(16) 



As we mentioned in the previous section, at the stage of 
oscillations when <j) M p the average value of the cur- 
vature R vanishes, so one can neglect the term ~ £cf> 2 R. 

The self-interaction \\4> A also leads to the generation 
of fluctuations of the field <j>. The equation for the eigen- 
modes 0k (i) is 



3- <pk 
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+ 3A0 2 U fc =0 
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Note that this equation is identical to equation ( |l6| ) with 
g 2 = 3A. Therefore, the study of the fluctuations cf>k 
in the jXcf) 4 model is a particular case of the general 
equation for fluctuations (|lq). 

The physical momentum, p = ^ly, in equation ( |T6| ) 
is redshifted in the same manner as the background field 
amplitude, <p(i) = ^£y. Therefore, the redshifting of mo- 
menta can be eliminated from the evolution of Xfc- In- 
deed, let us use the conformal transformation of the mode 



4 



function Xk(t) = a(t)xk(t) and rewrite the mode equa- 
tion for XjJt) with the dimensionless conformal time x 
(see Eq. (Ugh): 
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where for simplicity we drop the initial value of xq — 2.44. 
In this form the equation for fluctuations does not de- 
pend on the expansion of the universe and is completely 
reduced to the similar problem in Minkowski space-time. 
This is a special feature of the conformally invariant the- 
ory lA0 4 + i ff 2 2 X 2 . 

For the fluctuations of the field ip = acj> one has 
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Equation (|T|) will be the master equation for our in- 
vestigation of the resonance in the conformally invariant 
theory. The comoving momentum k enters the equation 
in the combination 



k 2 
\(p 2 



(20) 



Therefore the natural units of the momenta k is \f\(p. 
Equation ( fl8l ) describes oscillators, Xk, with a variable 
frequency 
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which periodically depends on time, x. It is well known 
that in this case the solutions Xk are exponentially 
unstable: Xk{x) oc e^ kX . If we choose the vacuum 
positive-frequency initial condition, X k (x) ~ e ,„ , we 



then expect the exponentially fast creation of x-particles 
(rife oc e 2fJ,kX ) as the inflaton field oscillates. The strength 
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of interaction with the periodic oscillations cn 2 (a 

is given by the dimensionless coupling parameter <? 2 /A. 
This means that the condition of a broad parametric res- 
onance does not require a large initial amplitude of the 
inflaton field, <po, as in the case of the quadratic po- 
tential Q. As we will see, the combination of param- 
eters g 2 /X ultimately defines the structure of the para- 
metric resonance in the theory. It turns out that the 
strength of the resonance depends rather non trivially 
(non-monotonically) on this parameter. 

From a mathematical point of view, the mode equa- 
tion (|l^) belongs to the class of Lame equations |12j . In 
the context of preheating this was first noticed in [Q, 
and then thoroughly studied for 0(N — ► oo) theory (i.e. 
for g 2 = A) in and ||. In this paper we perform a 
numerical and analytical investigation of the parametric 
amplification of fluctuations in the conformally invariant 

theory jXcf) 4 " + ^g 2 ip 2 x 2 for an arbitrary parameter ^r. In 
the next section, we present the two dimensional chart of 



the stability/instability bands for the Lame equation dig ) 

2 

in terms of variables k 2 and In subsequent sections, 
we give a new analytic treatment of the Lame equation 
in the case ^- = n (" +1 ) with integer n. We will also 
perform an analytical investigation of the resonance for 
X < 1 and for — > 1. 



IV. STABILITY/INSTABILITY CHART IN THE 
CONFORMAL THEORY 

As we shown in the previous section, the equation for 
vacuum fluctuations interacting with the inflaton oscil- 
lations in the conformal theories can be reduced to the 
similar problem in the Minkowski space. The equation 
for fluctuations (11^) in this case contains only two pa- 

2 

rameters. The first parameter is which which features 
the strength of the interaction. The second parameter is 
the momentum of vacuum fluctuations k in units of the 
frequency of the inflaton oscillations. As is well known, 
the solutions Xk of this equation may be stable or un- 

2 

stable depending on the particular values for k and 
considered. At the stage of the free resonance when we 
do not take into account the backreaction of the unsta- 
ble fluctuations, Eq. ( |l8| ) is an equation with periodic 
coefficients, which belongs to the class of the Lame equa- 
tions. The stability/instability chart of another equation 
with periodic coefficients, the Mathieu equation, is well 
known and can be found in many textbooks, see e.g. 
We are unaware of the stability/instability charts for the 
Lame equation, which describes preheating in the con- 
formally invariant theories. Therefore in this section we 
present the stability/instability chart Eq. ( |l8| ) in vari- 
ables (k 2 , 5r)' wruc h wc obtained by solving this equa- 
tion numerically. 

Fig. |^ shows a typical resonant solution of equa- 
tion (|l^). Though we have plotted the particular case 

k 2 = 1.6, ^- = 3, the form of the resonant solution is 
generic. The upper plot demonstrates the amplification 
of the real part of the eigenmode Xk {x) (solid curve) in 
the oscillating inflaton background (dotted curve). 
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FIG. 3. The typical resonant production of particles at the 
particular choice of rescaled comoving momentum k 2 = 1.6, 

2 

and the parameter 4- = 3. The upper plot shows the amplifi- 
cation of the real part of the eigenmode Xt{x). The lower plot 
shows the logarithm of the comoving particle number density, 
rife, calculated with formula (|2^). The number of particles 
grows exponentially, log rife ~ 2[ikX- In this case, Hk ~ 0.035. 

In addition to the investigation of the rapidly oscillat- 
ing functions X^{x), it is convenient for analytical and 
numerical work to consider the evolution of the comoving 
number density of created x-particles, rik, with comov- 
ing momentum k. This can be defined from the comoving 
energy density and the energy per particle, u>k'- 

The lower plot of Fig. | shows the evolution of the log- 
arithm of rjfc (solid curve) and the inflaton field (dotted 
curve). For the growing solutions after an initial tran- 
sitional period the number of particles increases expo- 
nentially, Inrifc 2/i fc x, where /X& is the characteristic 
exponent of the unstable solution. In the particular case 
shown, /ifc sa 0.035. 

2 

For arbitrary values of n and we can obtain a nu- 
merical solution of equation (|l8|) and exploit the simple 
relation mn fc ss 2/i fc a; to extract the characteristic expo- 
nent for the growing modes. For the regions of stability 
the characteristic exponent formally is imaginary. In this 
way, the stability/instability chart for the Lame equation, 
Fig. ||, is constructed. Shaded (unshaded) regions of the 

chart indicate values of k 2 and ^- for which the solutions 
are unstable (stable). For the instability bands, a darker 
shade indicates a larger characteristic exponent. An im- 

2 

mediate result is that, for a given range of the largest 
characteristic exponent will occur for k 2 = between the 
integer values ^- = "( n 2 +1 ) with n integer. 
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FIG. 4. The stability/instability chart for the Lame equa- 
tion for fluctuations Xk(x) in the variables (k , ^-), obtained 
from the numerical solution of equation (|18|) . Shaded (un- 
shaded) areas are regions of instability (stability). For insta- 
bility bands, the darker shade implies a larger characteristic 
exponent ^ife. Altogether, there are 10 color steps. One color 
step corresponds to the increment A/ife = 0.0237, so the dark- 
est shade corresponds to maximal fit = 0.237, the least dark 
shade in the instability bands corresponds to /ife = 0.009. For 
positive k 2 , there is only one instability band for the particu- 

2 

lar values of the parameter 4- = 1 and 3. This occurs because 

2 

the higher bands shrink into nodes as 4- approaches 1 and 3. 

This is demonstrated in Fig. ||, where slices of the sta- 
bility/instability chart show the characteristic exponent 

as a function of k 2 for various values of The top 

panel of Fig. | plots the cases ^ = 1.0, 1.5, 2.0, 2.5, 3.0, 

2 

labeled a through e respectively, v = 1 corresponds to 

2 

n = 1, t- = 3 corresponds to n = 2. As claimed, we 
see that the largest value of the characteristic exponent 
occurs for k 2 — at a value of ^- between the limits 1 
and 3 (curve c) . Similarly, the lower panel of Fig. |5| plots 
the cases si = 6.0, 7.0, 8.0, 9.0, 10.0, labeled a through e 

2 

respectively. The values 4- = 6 and 10 correspond to 
n = 3 and 4. Again we see that the largest value of the 
characteristic exponent occurs for k 2 = at a value of 

2 

3- between the limits 6 and 10 (curve c). 

This stability/instability chart is very similar to the 
stability/instability chart of the Mathieu equation, but 
there are important differences as well. For the Math- 
ieu equation there are infinitely many instability bands 
corresponding to each value of the parameter q, which is 

2 

analogous to our parameter ^- . Meanwhile for the Lame 
equation some of the instability bands may occasionally 

2 2 

shrink to a point. As a result, for ^- = 1 and for ^- = 3 
(Fig. ||, curves a and e respectively) there is only one 
instability band. This will be shown analytically in sec- 
tions VI and VII. From the stability/instability chart for 
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the Lame equation, Fig. Q, we see that, topologically, 
this occurs because the higher zones shrink to nodes as 

2 

4- approaches 1 and 3. 





FIG. 5. Slices of the stability/instability chart, Fig. [| , 
reveal the dependence of the characteristic exponent, [ik, on 
k for several particular values of v- For top panel <h- = 1.0, 
1.5, 2.0, 2.5, and 3.0, labeled a through e respectively. The 
numerical curves a and e for \- = 1 and % = 3 are identical 
to the analytic predictions 
Section VII. For lower panel - x 
labeled a through e respectively. 



of Section VI and © of 
= 6.0, 7.0, 8.0, 9.0, and 10, 



Similarly, there are a finite number of instability bands 

_ n(n+l) 
2 



for positive k 2 whenever C- = 



However, as for 

2 

the Mathieu equation, all other values of have an in- 
finite number of instability bands. This is true in par- 

2 

ticular for 4- <C 1 where, as we will show in Sec. VIII, 
the Lame equation may be formally transformed into the 
Mathieu equation with the parameters A w 1.3932k 2 and 
q « 0.3464^ < 1. 



and in the limit j « 1 



Thus, with this change of variables 

the stability /instability chart 
for the Lame equation is found to coincide exactly with 
that for the Mathieu equation. 

We now proceed to develop the new analytic results for 
preheating in the theory j<ft 4 , which as we have seen, are 
hinted at by the stability/instability chart for the Lame 



equation, Fig. |j. 



V. ANALYSIS OF THE EQUATION FOR 
FLUCTUATIONS 

In this section we begin the analytic investigation of 
the Lame equation ( |l8| ) for the fluctuations Xk{x). In 
particular, in the next two sections we will try to find 

2 

the values of the parameter 4- for which an analytical 
solutions can be obtain in closed form, and construct 
these solutions. 

We will also investigate the resonance in two limiting 

2 2 

cases: ^- <C 1 and 4- ^> 1. In the first case one can use 

2 

perturbation theory in the small parameter ^- < 1, see 

Section VIII Jn|. In the opposite limit, ^- 3> 1, we can 
implement the method of successive parabolic scattering 
Q , see Section IX. 

It is known that the Lame equation can be solved in 
terms of the transcendental Jacobi functions, which in 
turn are given by series expansions. Earlier we reported 
the result for the characteristic exponent /i = 0.0359 for 
A</> 4 theory jl4],(|] . Analytic investigation of the resonance 
using these transcendental functions gives the width of 
the unstable zone and the maximum of the characteristic 
exponent, /ij., in the physically interesting cases of the 

2 

0(N — > oo) theory = 1 in our convention) and the 

l\(j} 4 self-interacting theory = 3) 

However, calculations involving these transcendental 
functions are extremely tedious. Fortunately, it turns 
out that for 



g 2 n(n + 1) 
T = 2 



(23) 



with n an integer, one can obtain simple, closed-form 
solutions to the master equation (18). This includes in 
particular the most interesting cases g 2 = A and g 2 = 3A. 
To find the solutions of the fluctuation equation dlg| ) 

for ^- = " ( -" 2 +1 - ) , we will rewrite Eq. ( |l8[ ) in the so-called 
algebraic form. We will use the "time" variable z instead 
of x: 



z(x) 



cn 2 ( x, 



1 



d 
dx 



^2z(l-z 2 ) 



dz 



(24) 



Equation ( p3| ) for fluctuations becomes 



d 2 X 



dX k 



2z(l - z 2 )—^- + (1 - 3z 2 )— -1 + (n 2 + *—z)X k = 
dz z \ l 



dz 



(25) 



Omitting the lower index k for simplicity, let Xj(z) and 
Xi{z) be two linearly-independent solutions of (p5|). One 
of them exponentially grows, another exponentially de- 
creases during the resonance. Let us also introduce the 
bilinear combinations X 2 7 X 2 , and XiX2- From (|25|) 
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it follows that these bilinear combinations obey a third 
order equation 



„ , o ^d 3 M . 2 d 2 M 
2z{z 2 ~\) — - + (9z 2 -3)- 



dz 3 

2y -3)z + 2n 2 



dz 2 

dM q 2 

— 2^-M = 

dz A 



(26) 



The three solutions, M(z), of this equation correspond 
to the three bilinear combinations of X\ and X 2 . The 
crucial observation is that for 4- = the equation 

( p6| ) admits a polynomial solution of degree n. In the 
particular cases n = 1 and rc. = 2, we have 



ra=l: Mi(z) = z-2k 2 

.2_ 2 

3 



n = 2 



M 2 (z) = z 2 



-K 2 Z - 1 + -K 4 



(27) 



Obviously, the polynomial function M(z) must be the 
product of an exponentially growing solution and an ex- 
ponentially decreasing one, i.e. M(z) = Xt{z)X^{z) in 
the resonance zone. From this, as we will show in the 
next two sections, one can construct the closed-form so- 
lutions X(z). 

Therefore, in the physically interesting cases n = 1 and 
ii = 2 we will obtain simple closed form solutions instead 
of the complicated transcendental functions. This signif- 
icantly simplifies the study of preheating in these cases. 
In particular, we will find the form of the characteristic 
exponent [i^ as a function of k 2 in each case. 



VI. CLOSED FORM SOLUTION FOR 



1 



In the case g 2 = A equation ( pq ) in the resonance band 
gives 



X 1 {z)X 2 (z) = M x {z) , 

where 

M x (z) = z-2k 2 . 
The Wronskian of equation (pq) for X(z) is 

C 



dX 2 dX, 

Ai— A2 



dz 



dz y/ Z (l-Z 2 ) ' 



(28) 
(29) 

(30) 



where C is some constant, C = Ci, to be defined. From 
( p8| ) and (|30| ) we immediately obtain the closed form so- 
lutions 



X 1)2 (z) = v^MiMexp J -j= 



dz 



z 2 )M x (z) 



(31) 



Now, substituting this solution back into equation (|2 
for X(z), we find the constant C\\ 



Ci = ^2k 2 (1-4k 4 



(32) 



For exponentially growing solutions, Ci must be real; 
therefore the exponentially growing solutions for fluctu- 
ations with k 2 > take place in a single instability band 
for which 

< k 2 < I . (33) 



The growing solution of 



e fJ,kX P[z(x)], where P[z{x 



2 



has the form X(x) = 
a periodic function of the 
conformal time x. Using (|31j), we can now find the char- 
acteristic exponent fi^ as a function of k. The technical 
details can be found in the Appendix. 
The final answer is 



where an auxiliary function is 



t/2 



J(«) 



,1/2, 



1 + 2k 2 sin ( 



(34) 



(35) 



o 

Recall that T w 7.416. Eq. @ is one of the most im- 
portant analytic results of our paper. Some numerical 

2 

values of /i/. as function of «r for ^- = 1 calculated with 
(^J) are listed in the upper half of the table below. The 
analytic form ( |34| ) is in excellent agreement with the nu- 
merical results for this case plotted in the top panel of 
Fig. as curve a. The maximum value of the character- 
istic exponent for ^- = 1 is [i max ~ 0.1470 at n 2 fa 0.228, 
in agreement with the numerical value for [i max of Fig. 





K 2 




1 


0.0 


0.000 


1 


0.1 


0.1238 


1 


0.2 


0.1460 


1 


0.21 


0.1466 


1 


0.22 


0.1469 


1 


0.228 


0.1470 


1 


0.23 


0.1470 


1 


0.24 


0.1468 


1 


0.25 


0.1465 


1 


0.3 


0.1411 


1 


0.4 


0.1117 


1 


0.5 


0.000 


3 


1.5 


0.000 


3 


1.55 


0.02981 


3 


1.60 


0.03570 


3 


1.61 


0.03595 


3 


1.615 


0.03598 


3 


1.62 


0.03594 


3 


1.625 


0.03583 


3 


1.65 


0.03427 


3 


1.70 


0.02460 


3 


1.732 


0.00 
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VII. CLOSED FORM SOLUTION FOR 



= 3 



VIII. SOLUTION FOR % < 1 



The method of obtaining a closed-form analytic solu- 
tion, Xk(z), in the case g 2 = 3A is similar to that of the 
previous section. In the resonance zone with g 2 = 3A, 
equation (Eq) gives 



X 1 (z)X 2 (z) = M 2 (z) , 



where now 



M 2 {z) 



2 2 

— K Z - 

3 



1 + — K 

9 



(36) 



(37) 



The Wronskian of equation (|25| ) is the same as in ( |30[ ), 
but with a new constant, C = C 2 , to be defined for this 
case. Therefore, the closed form solutions are the same as 
in (|l|), but with M 2 (z) in place of Mi(z). Substituting 
this solution into equation (Eq), we find the constant C 2 
in this case 



'WiM"*-!)"-" 



(38) 



Therefore, in the case 



3 for k 2 > 0, there is also 



only a single instability band corresponding to 



< K < 



(39) 



For illustration, we plot the resonant solution Xj-(x) 
in the top panel of Fig. ||. Notice that Xk(x) oscillates 
twice within one inflaton oscillation. Using the solution 
( |3l| ) with M 2 (z) and C 2 , we can find fik in this case; see 
the Appendix for details. 



The resulting characteristic exponent for 4- = 3 is 



fJ-k 



8V2 
~9T~ 



^Jk 2 (k 4 



(3 



J(k) 



(40) 



where the auxiliary function J(k) is 



tt/2 



J(k) 



d6- 



|k 2 sin( 



+ (|k 4 - l)sin 2 6> 



(41) 



in this case. Formula (40) is another important result 
of our paper. Some numerical values of fj,k as a function 
of k 2 for \ = 3 calculated with @ are listed in the 
lower half of the last table. The analytic form (^) is in 
agreement with the numerical results for this case plotted 
in the top panel of Fig. [5] as curve e. The maximum 
value of the characteristic exponent for ^- = 3 is fi max ~ 
0.03598 at k 2 w 1.615, in agreement with the numerical 
value for fi max of Fig. |[ 



In this section we investigate the equation for fluctu- 
ations (|l^) in the limiting casejy- <C 1. Let us recall 
that f(x) is g iven by the series (ffq), and hence, f 2 {x) in 
equation ( [l8|) can be decomposed as 



f 2 {x) =F + F lC os( 



47ra;\ 



F 2 cos(— J 



(42) 



where F = 0.4570, F 1 = 0.4973, F 2 = 0.04290 and so 
on, but Fk = l- O ne can see k -X"*;^) in the form of 

a harmonic series of terms cos y ™ j with slowly varying 

2 

coefficients. If ^j- is a small parameter, one can develop 

2 

an iterative solution with respect to v- It is easy to 
show that the leading contribution to Xk{x) comes from 
the lower harmonic: cos(-^?). Keeping only this term, 
the equation for Xk{x) can be reduced to the Mathieu 
equation 



d 2 X k 
dr 2 



{A + 2qcos2r) X k = 



(43) 



where r= 2 f,A= (f£) 2 , and q = £ {^f F x . Thus, 
our theory is effectively reduced to the Mathieu equa- 
tion only in the limit q <C 1, where it has instabili- 



ties in very narrow resonant bands around k 2 = 2 -|? i , 

to = 1, 2, ... . The results of the numerical investigation 

of the instability zones plotted in Fig. indeed show that 

2 _ 1-1 

for ^- <C 1 the parametric resonance corresponds to that 

of the Mathieu equation. 




FIG. 6. The maximum value of the characteristic expo- 
nent Umax extracted from the stability/instability chart, Fig. 

^| is plotted as a function of ^- (solid curve). The func- 

2 

tion fM m ax(^) is non-monotonic. The universal upper limit 

of Umax is 0.2377. The local minima of the function are grad- 

2 

ually increasing with 4-, and asymptotically approach 0.2377. 

2 2 

The dotted line is the prediction fimax ~ 0.1467^- for y < 1, 
when the mode equation ( [l8| ) is effectively reduced to the 
Mathieu equation 
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The exponentially growing solution of the Mathieu 
equation, Xk(x) oc e^ kX , has a maximum characteristic 
exponent (in the first zone) 



El 

4A 



F 1 w 0.1467 



9 



(44) 



In Fig. p we plot the maximum value of the characteristic 

2 

exponent as a function of together with the prediction 
( |44| ) for n m ax from the Mathieu equation. As one can see 
from Fig. Eq. ( f44| ) works extremely well even up to 



IX. ANALYTIC SOLUTION FOR ^ > 1 

In this section we consider the limiting case when the 
parameter 4- is very large. In Fig. [7] we plot the time 
evolution of fluctuations X^(x) is this case. 



A', 




4 
3.5 

3 
2.5 

2 

lllTlfc 1.5 

I -. 

0.5 


-0.5 
-1 







FIG. 7. The same as in Fig. pbut for a large value ^- 3> 1, 

here for the particular choice °y- = 5050 and k 2 = 29.0. The 
upper plot shows the time-dependence of the real part of the 
eigenmode Xk(x), which demonstrates the adiabatic (semi- 
classical) behavior between zeros of the inflaton oscillations 
(dotted line), where the comoving occupation number nt of 
created particles is constant (lower plot). The lower plot 
shows log rik as a function of time x. Particle creation oc- 
curs in a step-like manner only in the vicinity of the zeros of 
the inflaton field, where the adiabaticity is broken. The en- 
velope of log rife is approximated by 2[ikX. The characteristic 
exponent for this example is fik ~ 0.1. 

In Fig. |7 we plot the number of particles rik[x) in 
a given mode as a function of time x calculated from 
Xfc(x) with Eq. (^). The basic observation is that, for 

2 

4- 3> 1, the evolution of the modes Xk(x) is adiabatic 
and the number of particles Uk(x) is constant between 
the zeros of the background field. Changes in the num- 
ber density of particles occur only near times x = Xj 
when the amplitude of the inflaton field crosses zero, i.e. 
tp(x — Xj) — 0. To describe the effect of a single kick at 
x = Xj, it is enough to consider the evolution of Xk(x) 
in the interval when <p 2 (x) is small and can therefore 
be represented by its quadratic part oc (x ~ Xj) 2 . This 
process looks like wave propagation in a parabolic poten- 
tial. Outside of these time intervals, X^(x) has a simple, 
semiclassical (adiabatic) form. We can combine the ac- 
tion of the subsequent parabolic potentials to find the 
net effect of particle creation. This method of succes- 
sive parabolic scattering was formulated and applied to 
the broad parametric resonance for the quadratic infla- 
ton potential in H. This method, as we see, can also be 

applied to the conformally invariant theory for 4- ^> 1. 

We expect that the semiclassical solution is valid ev- 
erywhere but around Xj . Thus, prior to scattering at Xj , 
the mode function Xf.(x) has the adiabatic form 



J UJk< ^ x 




(45) 



where the coefficients cxL and pL are constant for Xj-\ < 
x < xj and normalization yields \otk\ 2 — |/?fc| 2 = 1. After 
scattering when x = Xj, Xk(x) in the interval x-j < x < 
Xj+i again has the adiabatic form of equation (W3) but 
with new constant coefficients, cP^ 1 and ■ 

The form is essentially the asymptotic expression of 
the incoming waves (for x < Xj) and similarly for the the 
outgoing waves (for x > Xj) scattered from a parabolic 
potential (x — Xj) 2 at the moment Xj. Therefore, the 
outgoing amplitudes, c4 +1 and , can be expressed 
in terms of the incoming amplitudes, o? k and 0L, with 
the help of the reflection and transmission amplitudes for 
scattering at a parabolic potential g). For this we need 
the mode equation around a single parabolic potential at 



In the vicinity of Xj, cn 
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Then equation ( |18| ) around Xj is reduced to the simple 
equation 



d 2 X k ^ ( 2 u 
dx 2 



(46) 



The mapping of a k , [3i into cP k +1 , m terms of 

parameters in equation (Era) reads 



J+ 1 



vT 



e _ " e e 



le 2 




'\/2Cp 2 g 



where Cfc = argT ( 1+ * e and e 2 

The phase accumulated by the moment Xj is 9 k — 

T/i i 

J dx u)k{x) = j 9k, where 6k = 2 J dx J n 2 + \f(x) 
o o 
is the phase accumulating within half of a period of the 
inflaton oscillation. 

In the regime when a large number of particles have 
been created, n J k = \fi 3 k \ 2 ^ 1' we have \a J k \ « |/3£|, so aj, 
and (3 k are distinguished by their phases only. In this case 
there is a simple solution of the matrix equation (J47|) : 



V2 



V2 



(48) 



where •& is a constant phase and /!& is the characteristic 
exponent. T sa 7.416 is the period of oscillations of the 
inflaton field in the variable x, so the number of particles 
grows as e 2Mfc2 \ Another solution is similar to dlgj ) but 
with the substitution 9k — > Ok + t- 

Substituting the solution ( jig] ) into equation (fi7|), we 
get an equation for the parameters jik and 9k 



= |cos(6» fe -Cfc)|v / l + e-' re2 

+ ^(l + e-^^cos^efe-CO-l 



(49) 



In the instability zones, the parameter fif. of equation 
(fl9|) should be real. From this we obtain the condition 



tan( 



Cfe)l < e- 



(50) 



for the momentum k to be in a resonance band. 

To further analyze the conditions for the strength j50| ) 

and widths ([49] ) of the resonance, one should calculate 

2 

the phase 9k — Cfc- For ^- > 1 we have 



Go 




^/ 2 (x)-argr(i±^ 



(51) 



Using equations (|5l|), (p0|), and (|49|), we find the char- 
acteristics of the resonance in the regime ^> 1. From 
( |50| ) it follows that the resonance is efficient for e 2 < 7r _1 , 
i.e for 



K < 



2ir 2 \ 



(52) 



Equation (|5C|) transparently shows that, for a given 
there will be a sequence of stability/instability bands as 
a function of k. The width of an instability band, where 

the resonance occurs, is Ak 2 ~ v/fx- Let the integer 

be Z. From (BlJ) it follows 



part of the large number 

that if we vary n 2 within the range 2-7T 2 y^3j(ln ^r) *> 
then within this interval of k 2 the phase 9k — Cfc reaches 
either lir or (Z + l)7r. Then within this resonance band 
we get the maximum value [i m ax defined by the equation 
© with |cos(0 fe -Cfe)l = 1: 



e 2 ' 



+e"— . (53) 

The characteristic exponent /U TOa:c is a non-monotonic 

2 2 
function of If the value of the parameter 4- is ex- 
actly equal to 2Z 2 where Z is an integer, then the strongest 
resonance occurs at k 2 — 0,[] and from ( |5"3| ) we get 



// m ax = yln(l + V2) w 0.2377. 



(54) 



This is actually a general result for the upper limit of 



Umax for an arbitrary see Fig. [| 



If is not 



exactly equal to 2Z , then [i max occurs at a non-zero 

k 2 and is smaller than 0.2377. It is interesting that 

2 / 2 \ 

in the formal limit ^- — > 00 the function /i max (^-J 

asymptotically approaches the value 0.2377 for arbitrary 

2 

4-. To see this, we have to check that a variation of 



n" ~ \J 2n' 2 \ ( m 5r) 1 ^ S compatible with the condi- 
tion for an efficient resonance, e 2 < ir^ 1 . This occurs 
for ^ > e 27r2 

mal value of as a function of ^- very slowly increases 

towards 0.2377. Therefore, although fj, max is not a mono- 

2 2 

tonic function of 5r , for 5r S> 1 the resonance is stronger 



10 8 . In Fig. ^| we see that the mini- 



X ' A 

both in terms of the characteristic exponent /i„ 
width k 2 . 



and the 



X. BACKREACTION OF CREATED PARTICLES 

Thus far we have considered the parametric resonance 
in the conformally invariant theory (hi) in an expanding 



*It is easy to see from @ that the mode k 2 = is within 
the resonance band if 2l 2 



I < fi_ < 2Z 2 +1. 
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universe neglecting the backreaction of the amplified fluc- 
tuations of the fields and x.- I n the next two sections 
we will study the effects related to the backreaction. 

In the theory |A</> 4 + \g 2 <t> 2 X 2 : the equation of motion 
for the inflaton field <j>(t) looks as follows: 

4> + 3H 4> + X(f> 3 + Z\{<t> 2 )(t> + g 2 {x 2 )4> = . (55) 

The two additional terms are due to the one-loop Hartree 
diagrams; (c/f> 2 ) and (x 2 ) stand for the quantum fluc- 
tuations of the fields <j> and x respectively. There are 
also higher-loop corrections such as 2g 2 ((f>x 2 ) and 6A(</> 3 ), 
which are not necessarily negligible at the end of preheat- 
ing when we may expect n x ~ I/5 2 and ~ 1/A (or 
even higher, if one takes into account rescattering) . Here 
we will work in the one-loop approximation. 

If again we use the conformal transformation 77 = J — , 
Xk = and similarly <pk = <*Pk/o>i then the equation for 
the background field ipirf) is 

^" + \>f 3 + 3\((p 2 )(p + g 2 (X 2 )<p = , (56) 

with the comoving vacuum expectation values for X and 
ip correspondingly 

(57) 

The integral of the equation (^6|) coincides with the 
energy density 

Ptot = \<p' 2 + + P v + px ■ (58) 

The first two terms describe the energy of the classical 
field ip, p v and px correspond to the energy density of 
(^-particles and X-particles respectively: 

P v = -Jtz J d 3 k V^ 2 + 3A^ 2 nl , (59) 

Px = 7A3 / <? k Vk 2 + g 2 <P 2 nf . (60) 

Here nt and n£ correspond to the occupation numbers 
of the (/3-particles and X-particles. It is easy to show that 
p' = g 2 {x 2 )f l f' an d p' v = i\{(p 2 )(fip' ', and therefore Eq. 
(p§|) is an intergal of Eq. ([56"|). 

To close the set of self-consistent equations we need 
the equations for the modes <pk(%) and X/.: 

Vk(v) + (k 2 +Il v + 3\v 2 ( v ))^k = 0, (61) 

X'kiv) + {k 2 +Tl x +g 2 i p 2 (r ] ))Xk = . (62) 

The polarization operator H v consists of = 3\(ip 2 ) + 
g 2 (X 2 ) and the non-local term n 2 which emerges in the 



one-loop approximation beyond the Hartree diagram, see 
Fig. I 



X ip 




<p(t) <p(t) (p(t) (p(t) 



FIG. 8. The diagrams for the polarization operator of the 
field ifk- Thin and thick lines represent the fields 4> an d X 
respectively. Vertical lines correspond to the oscillating back- 
ground field <p(t). Tl v corresponds to the Hartree approxima- 
tion which takes into account the contribution of (X 2 ) and 
{<p 2 }- The contributions of II J, and 11^ to the effective mass 
of (^-particles can be comparable to each other. 

The calculation of the polarization operator n 2 in the 
regime of parametric resonance is rather complicated. 
Estimates of n 2 performed in ||] indicate that it can be 
of the same order of magnitude as the standard Hartree 
polarization operator n*. The polarization operator n 2 
was not taken into account in the previous treatment of 
the self-consistent equations for the eigenmodes in the 

approximation , but in fact it may survive in the 
limit N 00 j2|. 

Similarly, the polarization operator Hx is equal to 
g 2 (ip 2 ), plus an additional non-local term H x - We ex- 
pect that n v > 0, > 0, as suggested by the Hartree 
approximation. 

A complete calculation of the polarization operators 
n v and Hx is outside the scope of this paper. Fortu- 
nately, as we will see in the next section, one need not 
really know exact expressions for Tl v and Hx in order to 
make an estimate of the density of produced particles at 
the time when the feedback of the amplified fluctuations 
terminates the parametric resonance. 

XI. DYNAMICAL RESTRUCTURING OF THE 
RESONANCE 

In this paper we found that the structure of the para- 
metric resonance in terms of its strength and width 
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strongly depends on the parameters of the model. For 
example, the parametric resonance in the simplest con- 
formally invariant theory jXff) 4 +\g 2 cj) 2 x 2 is very different 



from that in the theory ^m^cj) 2 



In the sim- 



_ 

plest conformally invariant theories which we consider in 
this paper the structure of the resonance is determined 
by the combination g 2 /X. 

How does the resonance develops if the backreaction 
of the accumulating fluctuations is taken into account? 
The answer to this question also strongly depends on the 
parameter g 2 /X. 

For illustration we consider the model of the self- 
interacting inflaton field \X(j) A 1 no \ field is involved. In 
this case we shall take g 2 = in all the equations (56), 
( |6l"| ) , d62| ) . As we already mentioned, if one neglects back- 
reaction, the equations describing the resonance for the 
modes tpk in this theory coincide with the equations for 
the modes Xk in the theory with g 2 = 3A. Thus, we 
can use the results of the investigation of the theory with 
g 2 = 3A obtained in Sec. VII for our analysis. 

Historically, the model jX<p 4 was one of the first mod- 
els illustrating the general idea of preheating. The in- 
vestigation of the stability/instability chart for the Lame 
equation has shown that this model in a certain sense is 
the least favorable for the development of the resonance: 
it has only one resonance band, and the characteristic ex- 
ponent fi for this theory is anomalously small, see Figs. 
4, 6. Originally it was expected that preheating in this 
model would rapidly transfer about half of the energy of 
the oscillating scalar field to the (^-particles, after which 
the decay of the field <fi would continues at a much slower 
pace. However, the results of computer simulations of 
preheating in this theory indicated that the stage of effi- 
cient preheating ends as soon as the fluctuations of pro- 
duced particles (if 2 ) grow to 0.05 <p 2 |5j. The interpre- 
tation of this result, however, was not quite clear. It 
was conjectured that the resonance terminates because 
of rescattering of the 0-particles. It was not clear also 
whether the decay of the field 4> continues at a slower 
pace until this field completely decays, or its decay even- 
tually shuts down. 

A complete investigation of this issue is rather difficult. 
First of all, the theory of rescattering is not fully devel- 
oped: various approximations often break down near the 
end of preheating when the occupation numbers of parti- 
cles are anomalously large (nk ~ A -1 ) 0,0. Even in the 
Hartree approximation (or in the 1 /TV-approximation) an 
investigation is extremely complicated |QH because it is 
very difficult to work with the solutions of equations for 
the growing modes in terms of the transcendental Jacobi 
functions. It may be easier to work with the solutions 
obtained in Sec. VI and VII. We will not perform a 
complete investigation of this issue here because, as we 
argued in the previous section, one may need to calculate 
the polarization operator beyond the Hartree approxima- 
tion, see |Q. Instead, we will make some simple estimates 
which will allow us to elucidate the mechanism which ter- 



minates the resonance in the theory \X<j) 4 . 

As we will see, the main reason for the termination of 
the resonance in the theory |A</> 4 is the restructuring of 
the resonance band due to the backreaction of created 
particles. This process occurs at (ip 2 ) <C ip 2 because 
the resonance band is very narrow. In the beginning 
of preheating in the theory jXcfi 4 the instability band is 
given by the condition 1.5A^5q < k 2 < 1.73A</5q, where (po 
is the initial amplitude of the oscillations of the field ip 
(|39| ) . It is sufficient to shift the position of the resonance 
band in momentum space by few percent, and the leading 
resonant modes Xk which have been growing since the 
start of the parametric resonance will not grow anymore. 
This will effectively shut down the resonance. 

There are two different effects which lead to a restruc- 
turing of the resonance band, and these effects act in 
opposite directions. First of all, particle production re- 
duces the energy of the scalar field, and therefore reduces 
the amplitude of its oscillations. This effect tends to re- 
duce the frequency of the oscillations and to move the 
resonance band towards smaller k. On the other hand, 
the effective mass of the field ip grows due to its inter- 
action with the (^-particles. This effect increases the fre- 
quency of oscillations and tends to shift the resonance 
band towards larger k. We will consider here both of 
these effects. 

To investigate the decrease of the amplitude of the os- 
cillations due to particle production, one should compare 
the total energy of the system before and after the ap- 
pearance of (</> 2 ): 



A -4 
7^o 



A 



1 



(27T) 



d z k^/k 2 + 3Xp 2 n k . (63) 



Here we calculate the energy density at the moment when 
ip' = 0, and the oscillating field is equal to its amplitude 
<p(r]). This amplitude is smaller than ipo due to the trans- 
fer of energy to the created particles. 

The resonance is most efficient in a small vicinity of 
k 2 1.6A<y5 2 . Therefore, the leading contribution to 
is given by integration near k 2 — 1.6A</3 2 : 



1 



(2tt) 3 



d 3 fcv/4.6A^ 2 n k = V4~6X<pn,p . (64) 



Eqs. (B3) and (64) give 



tp(rj) « ip 



V46An^ 



V 2 



(65) 



Thus, the creation of ^-particles diminishes the fre- 
quency of oscillations, because the frequency of oscilla- 
tions of the field ip in the theory Xcj) 4 is proportional to its 
amplitude. To evaluate the significance of this effect one 
may express it in terms of (4> 2 ) calculated at <p(rj) = tp: 



1 



d 3 k 



(2?r) 3 J ^/k 2 + 3Xp 2 y/HXip 



(66) 
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From the last two equations one obtains 



1 - 9.2^ 



(67) 



which leads to a proportional shift of the resonance band 
towards smaller k 2 . This indicates that even a very small 
amount of fluctuations (4> 2 ) ~ 10 _2 </?g may shift the res- 
onance band away from its original position, which may 
terminate the resonance for the leading modes (fk- 

This effect is partially compensated by the growth of 
the effective mass of the field ip. We will analyse this 
effect in the Hartree approximation, in which the field ip 
acquires the effective mass squared = 3X(p 2 ). One 
may relate II^ = 3\((p 2 ) to the number density of tfi- 
particles in the following way: 



3A 



d 3 k n k (r]) 



3A 



d 3 k rik 



3\p 2 ( V ) 



(2tt) 3 J ^1.6^ 2 + 3A(^ 2 y/1.6<p 2 + 3\<p 2 (r 1 ) 



(68) 



Note, that this quantity is time-dependent. It oscil- 
lates; its magnitude changes considerably several times 
within a single oscillation of the inflaton field, and it also 
grows exponentially during the resonance. The number 
density of tp-particles also oscillates and grows exponen- 
tially, but typically its oscillations are less wild than the 
oscillations of (<p 2 )- In the first approximation, we will 
neglect the oscillations of n^{j]). Also, we are trying to 
find the time when the resonance terminates, and at that 
time the average number density of particles becomes 
nearly constant. It is still difficult to find an analytic 
solution for ip^, with the time-dependent polarization op- 
erator (|6^), but one can easily find the solution numeri- 
cally. 

The result of the combined investigation of the two 
effects discussed above shows that the resonance on the 
leading modes ipk effectively terminates as soon as (<fi 2 ) 
grows up to 



0.05, 



(69) 



Note that even after this moment the resonance may con- 
tinue for a while for the new modes which can be ampli- 
fied in the restructured resonance band. However, this 
process is much less efficient. Thus, in the pure A</> 4 the- 
ory the rapid development of the resonance ends when 
the dispersion of amplified fluctuations is about 20% of 
the amplitude of the inflaton field, which corresponds to 
only 0.2% of the total energy. This result is based on 
rather rough estimates neglecting rescattering. It is in- 
teresting, however, that it is in complete agreement with 
the result of the lattice simulation of the parametric res- 
onance in the theory A0 4 || . 

We should emphasize that there are several specific 
reasons why the resonance in the particular case of the 



theory A</> 4 is relatively inefficient. First of all, the reso- 
nance band in this theory is very narrow and the charac- 
teristic exponent ijl is extremely small. This is no longer 
the case when one considers, for example, the theory de- 
scribing a x-field with g 2 = A or with g 2 = 2A. In these 
theories the characteristic exponent is much greater, the 
resonance band is rather broad, and it begins at k = 0. 
As a result, it is much more difficult to shut down the 
resonance in such theories. 

In the theories with a massive inflaton field there is an 
additional effect which makes the resonance more stable. 
Broad parametric resonance in such theories is stochastic, 
which makes it more difficult to shut down j2|. Now 
we are going to study what happens to the resonance 
in the conformally invariant theories if this invariance is 
broken by a small mass term. As we will see, stochastic 
resonance may appear in such theories as well. 



XII. PREHEATING IN THE THEORY OF A 
MASSIVE SELF-INTERACTING INFLATON 
FIELD 



In our previous paper pi we investigated paramet- 

ric resonance in the theory ^-^> 2 + ^-0 2 X 2 - We have 
found that reheating can be efficient in this theory only 
if g$ ^> m, where $ is the amplitude of oscillations of the 
inflaton field. This amplitude is extremely large immedi- 
ately after inflation, <f> ~ 10 
as 



M p , and later it decreases 



M p 
3mt 



(70) 



Due to this decrease, the ratio — rapidly changes. As 
a result, the broad parametric resonance regime in this 
theory is a stochastic process, which we called stochastic 
resonance. 

Here we studied the theory j(f> 4 + ^(j^x 2 for various 
relations between the coupling constants g 2 and A. In this 
theory the amplitude of the field <j> also decreases in an 
expanding universe, but it does not make the resonance 
stochastic because all parameters of the resonance scale 
in the same way as <f> due to the conformal invariance. 
One may wonder, what is the relation between these two 
theories? Indeed, neither of these two theories is com- 
pletely general. In the theory of the massive scalar field 
one may expect terms ~ j<fi 4 to appear because of ra- 
diative corrections. On the other hand, in many realistic 
theories the effective potential is quadratic with respect 
to <j) near the minimum of the effective potential. 

To address this question, let us study the theory 

^j-(j) 2 + j<p 4 + ^4> 2 x 2 - One may expect that for <j> > ^ 
parametric resonance in this theory occurs in the same 
way as in the model j </> 4 + ^^X 2 ■> whereas for <j> -C 

the resonance develops as in the theory ^(f> 2 + ^-</> 2 X 2 - 
Let us check whether this is really the case, ignoring for 
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simplicity the effects of backreaction of created particles, 
which is always possible in the beginning of the resonance 
regime. 

First of all, one should remember that at the beginning 
of the stage of oscillations in this theory one has $ ~ 
10 _1 Af p . Therefore there are two basic possibilities. If 

71 » 10 



_1 M p , then the term j(f) A never plays any role 



in determining the frequency of oscillations of the field 4>. 
Also, in this regime the particles <f> are not produced by 
parametric resonance, because the condition y/X(j) > m 
(analogous to the condition g<p > m for the production 
of x-particlcs 0) is violated. In such a case x-particles 
can be produced if lQ~]gM p ^> to. The theory of this 
process is described in H; we do not have anything new 
to add here. 

Another possibility, which we are going to study here 
in more detail, is that <C 10 _1 M P . Then in the begin- 

2 9 

ning the mass term ^-(f> does not affect the frequency 
of the oscillating scalar field <f>. Therefore, one could ex- 
pect that as the amplitude <& decreases from 10 _1 M p to 
^7=, the theory of parametric resonance coincides with 
the one described in this paper. 

2 

However, for large t- the situation is more compli- 
cated. Even though the mass term for 10~ 1 M p > $ ^> 
^7= does not affect the frequency of oscillations, it may 
affect the nature of the broad parametric resonance by 
inducing an additional rotation of the phase of the modes 
Xk- 

The reason why the broad resonance in the theory 
I |-0 2 -|-^-0 2 X 2 was stochastic can be explained as follows. 
The %-particles are produced when the field <p(t) comes 
close to the point = 0, which happens once during 
each time period At = — . During this time the phase of 
each mode Xfe grows approximately by <7$(t)7rm _1 . Dur- 
ing the next half of a period of an oscillation it changes 
by g$(t+ ^irm- 1 sa g^^wm' 1 + g$(^7r 2 m~ 2 . This 
destroys the phase coherence required for the ordi- 
nary resonance and makes the resonance stochastic if 
\g^>(t)TT 2 m- 2 \ > 1. 

The condition for the stochastic resonance in the the- 
ory ^ 2 - 



4> X can De obtained from Eq. (yu 



$ > 



> 



(71) 



Note 



In particular, for $ = ^= it gives -j= 
that by our assumption 3> 1. 

2 r, 

The generalization of this result for the theory ^-(j) + 

j4> 4 + ^-0 2 x 2 is straightforward, but the result is some- 
what unexpected. As a rough estimate of the time At one 
can take tt(2A$ 2 + to 2 )- 1 / 2 = n(2\ip 2 a~ 2 (t) + to 2 )" 1 / 2 , 
where ip = $a _1 (t) is the time-independent amplitude. 
The phase shift during this time is given by gipir(2\ip 2 + 
TO 2 a 2 (t)) -1 ' 2 . Thus, for to = this quantity is time- 



independent, and one can have a regular stable reso- 
nance. In the limit $ ^> -^7= one can represent the phase 

shift as 77^(1 — m i\Jf > ) ■ The change in this shift during 
one oscillation is "f^" , where H = § = . This 

gives the following condition for stochastic resonance: 



9 2 

. fit 7T Til 

$ < _£ 

~ 3AM P ' 



(72) 



% it gives £ > ^ 



Again, for $ 

This conclusion is illustrated by Fig. 9, where we show 

the development of the resonance both for the massless 

2 

theory with 4- ~ 1700, and for the theory with a small 
mass to. As we see, in the purely massless theory the log- 
arithm of the number density nk for the leading growing 
mode increases linearly in time x, whereas in the pres- 
ence of a mass to, which we took to be much smaller 
than y/X(f> during the whole process, the resonance be- 
comes stochastic. 




FIG. 

4.2 



9. 



Development of the resonance in the theory 
^-<f>" + ^<j>* + ^(t> 2 X 2 for sL = 5200. The upper curve cor- 
responds to the massless theory, the lower curve describes 
stochastic resonance with a theory with a mass m which is 
chosen to be much smaller than ^f\<j> during the whole period 
of calculations. Nevertheless, the presence of a small mass 
term completely changes the development of the resonance. 

In fact, the development of the resonance is rather com- 

2 

plicated even for smaller ^- . The resonance for a massive 
field with to <§; \fXcf> in this case is not stochastic, but it 
may consist of stages of regular resonance separated by 
the stages without any resonance, see Fig. 10. 
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FIG. 
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10. 



Development of the resonance in the theory 



2-r +h 



+ V 



with m 2 « A<£ 2 for 2- = 240. In 



this particular case the resonance is not stochastic. As time 
x grows, the relative contribution of the mass term to the 
equation describing the resonance also grows. This shifts the 
mode from one instability band to another. 



Thus we see that the presence of the mass term 2=j- (f> 2 
can modify the nature of the resonance even if this term 
is much smaller than j(f> 4 . This is a rather unexpected 
conclusion, which is an additional manifestation of the 
nonperturbative nature of preheating. This subject de- 
serves separate investigation 

Different regimes of parametric resonance in the the- 
ory ^(f) 2 + j<j) 4 + — </> 2 x 2 are shown in Fig. 11. We 
suppose that immediately after inflation the amplitude 
<f> of the oscillating inflaton field is greater than If 
_g_ < Vam p 



the x-particles are produced in the regu- 



lar stable resonance regime until the amplitude $(t) de- 
creases to -^t= , after which the resonance occurs as in the 

theory ^-(j) 
stochastic. 

If 4= > 



\<t> 2 x 2 



[§. The resonance never becomes 



the resonance originally develops as 
4 ^ al^ x ^ but 



in the conformally invariant theory j(j> 4 + %> 
with a decrease of the resonance becomes stochas 



tic. Again, for $(t) 



< m 



the resonance occurs as in the 



theory y^' + v^X • I n au cases the resonance eventu- 
ally disappears when the field <&(i) becomes sufficiently 
small. As we already mentioned in , reheating in this 
class of models can be complete only if there is a sym- 
metry breaking in the theory, i.e. m 2 < 0, or if one adds 
interaction of the field <p with fermions. In both cases 
the last stages of reheating are described by perturba- 
tion theory |l6|JT^] . 



stable resonance 



stochastic 
resonance 



stochastic resonance 
stable resonance 



no 








FIG. 11. Schematic representation of different regimes 
which are possible in the theory ^4> 2 + j4> 4 + ^^X 2 f° r 
< 10" 1 M P and for various relations between g 2 and A 
in an expanding universe. The theory developed in this pa- 
per describes the resonance in the white area above the line 
<& = ^= . The theory of preheating for $ < is given in Q . 
A complete decay of the inflaton is possible only if additional 
interactions are present in the theory which allow one inflaton 
particle to decay to several other particles, for example, an 
interaction with fermions ipip(f>. 

Adding fermions does not alter the description of the 
stage of parametric resonance. Meanwhile the change of 
sign of m 2 does lead to substantial changes in the theory 
of preheating, see Fig. O. We will investigate preheating 
in the theory —^(f) 2 + j(f> 4 + — <^ 2 x 2 in a separate pub- 
lication jl8). Here we will briefly describe the structure 
of the resonance for various g 2 and A neglecting effects of 
backreaction. This will give us a more general perspec- 
tive on the theory of reheating. 

First of all, at $ ^ m/\/~\ the field (j) oscillates in 

the same way as in the massless theory j<f> 4 + ^<j) 2 x 2 - 
Moreover, the condition for the resonance to be stochas- 
tic remains the same as before: $ < -2= , see Eq. 



However, as soon as the amplitude $ drops down to 
^7=, the situation changes dramatically. First of all, de- 
pending on the values of parameters the field rolls to one 
of the minima of its effective potential at <f> = The 
description of this process is rather complicated. De- 
pending on the values of parameters and on the relation 
between y/ (tf> 2 ), \J (x 2 ) and a = -^=, the universe may 
become divided into domains with 4> = ±cr, or it may 
end up in a single state with a definite sign of (f>. We 
will describe this bifurcation period in After this 

transitional period the field <f> oscillates near the min- 
imum of the effective potential at 4> = ±^ with an 
amplitude $ <C a = . These oscillations lead to para- 
metric resonance with x-particle production which can 
be (approximately) described as a narrow resonance in 
the first instability band of the Mathieu equation with 



1G 



A k = 4 fe2+g a g2 , g = For definiteness we will con- 

sider here the regime A 3//2 M p < to <§; \}/ 2 M p . The reSO- 
rt 2 1 

nance in this instability band is possible only if < 
the resonance in higher instability bands is very ineffi- 
cient and rapidly shuts down due to the expansion of the 
universe. Using the results of pi one can show that the 

\ 2 

resonance in the first band also terminates at $ < /?. . 

g M p 

By taking the upper limit of this inequality at $ ~ ^= 
one concludes that this resonance is possible only for 
~/a ~ ( \/Xm ) 1 ^ 4, (The resonance may terminate some- 
what earlier if the particles produced by the parametric 
resonance give a considerable contribution to the energy 
density of the universe.) However, this is not the end 
of reheating, because the perturbative decay of the in- 
flaton field remains possible. It occurs with the decay 

4 

rate T(<j) — > \X) = f^x 1 - This is the process which is 
responsible for the last stages of the decay of the inflaton 
field. It occurs only if one ^-particle can decay into two 
%-particles, which implies that < ~. 




FIG. 12. Schematic representation of different regimes 
which are possible in the theory —^<f) 2 + j(f> 4 + 4r0 2 x 2 - 
White regions correspond to the regime of a regular stable 
resonance, a small dark region in the left corner near the 
origin corresponds to the perturbative decay <j> \\. Unless 
additional interactions are included (see the previous Figure), 
a complete decay of the inflaton field is possible only in this 
small area. 



XIII. DISCUSSION 

In this paper we investigated the development of para- 
metric resonance in the conformally invariant theories of 
the type of j<j) 4 + ^-</> 2 x 2 - We have found that the de- 
velopment of parametric resonance in these theories does 
not depend on the expansion of the universe, and can be 
classified in terms of the ratio g 2 IX. This ratio deter- 
mines the structure of the stability/instability bands for 
the equations describing the resonance. 



We have found that the behavior of the resonance with 
respect to x-particle production is a non-monotonic func- 
tion of g 2 /X. For example, for g 2 — X and for g 2 = 3A 
equation for the perturbations of the field x has only one 
instability band, for ^- = "^" 2 +1 ^ there is only a finite 
number of instability bands, whereas for all other values 

2 

of ^- the number of instability bands is infinite. 

It is interesting that ^-particle production is least ef- 
ficient for g 2 <C A and for g 2 = 3A. For example, the 
characteristic exponent /i max for g 2 = 2A and for g 2 = 8A 
is almost 7 times greater than /z max for g 2 — 3A, see Fig. 
6. Meanwhile the characteristic exponent for the pro- 
duction of ^-particles in the theory j(f) 4 coincides with 
that of the field x f° r <? 2 = 3 A. Therefore x-particle pro- 
duction is typically more efficient than the production of 
^-particles (unless g 2 <C A). The non-monotonic depen- 

2 

dence of \i on the ratio ^- suggests that there exists an 
"unnatural selection" rule: The particles which are espe- 
cially intensively produced during preheating are not the 
ones which have the strongest coupling to the inflaton 
field, but those for which the characteristic exponent fj, 
is the greatest. 

In the conformally invariant theories the expansion of 
the universe does not hamper the resonance, so it ends 
only due to the backreaction of the produced particles. 
There are several different mechanisms which may ter- 
minate the parametric resonance. First of all, creation 
of particles leads to a decrease in the amplitude of os- 
cillations of the field <p = a<f>, which otherwise would 
remain constant. This leads to a proportional decrease 
in the frequency of oscillations in terms of the conformal 
time 77, which may shift the position of the instability 
band towards smaller momenta. There is also an oppos- 
ing effect which increases the frequency of oscillations 
due to the interaction of the homogeneous inflaton field 
with the produced particles. Finally, quantum fluctua- 
tions of the fields <fi and \ acquire contributions to their 
masses, which changes their spectra. A combination of 
all these effects leads to restructuring of the instability 
bands. This terminates the amplification of the leading 
modes which have been growing from the very beginning 
of preheating. Additionally, one may envisage effects re- 
lated to rescattering of produced particles, which may 
terminate the resonance even somewhat earlier. In this 
respect it is interesting that our estimates ignoring the 
process of rescattering give results which are in a very 
good numerical agreement with the results of computer 
simulations of reheating in the theory X<\> 4 performed in 
H where all of these effects including rescattering have 
been taken into account. 

Rescattering may be more important for g 2 3> A ||[l0) • 
However, in this regime one may need to take into ac- 
count possible small mass terms which should be present 
in realistic versions of the theory. As we have found, for 
g 2 ^S> X these mass terms lead to a radical change in 
the structure of the resonance not at $ < m/VX, as one 
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could naively expect, but much earlier, at <I> < jjjj-- 
In this regime the resonance becomes stochastic, the ef- 
fective width of the resonance band increases, making it 
much more stable with respect to various backreaction 
effects including rescattering [||. 

We should emphasize again that preheating is but the 
first stage of reheating, which does not lead to a com- 
plete decay of the inflaton field in any models which we 
studied so far. The last stages of preheating are always 
described by the perturbation theory |lqj , which will be 
developed further in our subsequent publication [jlTj . To 
illustrate this point, we described the development of the 
parametric resonance in the general class of models with 

the effective potential V(<j>, x) = ± ,J j-(/) 2 + + ^4> 2 X 2 - 
We have found that in these theories (without any other 
fields being added) the inflaton field can completely de- 

cay only if the sign of the term ^<j) is negative, which 
corresponds to spontaneous symmetry breaking. More- 
over j_this process completes only for ^ <\, see Figs. [lT| 



andfl| 



A complete inflaton decay is possible for 4- > ^ as well, 
even without spontaneous symmetry breaking, but only 
if the inflaton field has some other interactions, such as 
an interaction with fermions {jjipcj) with mass m^p < m/\/2 
jlTj . This conclusion implies that the decay of the infla- 
ton field is by no means automatic even if it is heavy and 
strongly interacts with other fields. Generically, the infla- 
ton field accumulates an enormously large energy density, 
which can be completely released only if it interacts with 
other particles in a very specific way |l]]. 

To understand how these results may change our point 
of view on the thermal history of the universe, let us 
suppose for a moment that the inflaton field does not 
have any interactions with light fermions, and that it 
has an effective potential — ,j y4> 2 + j <^ 4 + \^X 2 with 
\ ~ g 2 ~ 10~ 13 and with a small mass m ~ 10 2 GeV 
protected by supersymmetry. Then the final stage of 
reheating of the universe will begin only after the sym- 
metry breaking in this theory, and the reheating tem- 
perature estimated in accordance with Q will be smaller 
than 10 2 — 10 3 GeV. In such a theory the electroweak 
phase transition may never happen, or it may occur in 
an entirely different way. From the end of inflation until 
the symmetry breaking and the final stage of reheating, 
the universe will remain far away from thermal equilib- 
rium, and various nonthermal phase transitions and ex- 
plosive processes of particle production may occur. In 
such a model one should reconsider all issues related to 
the primordial gravitino problem, moduli field problem, 
baryogenesis, etc. 

The main conclusion of our investigation can be for- 
mulated as follows. The first stages of the process of the 
inflaton decay may occur much more efficiently than was 
previously thought, due to the effect of parametric reso- 
nance. The last stage of this process may be completely 
inefficient even if the coupling of the inflaton field to mat- 



ter is very strong, or it may be efficient only in a very 
narrow range of parameters, see Figs, [n] and [l2|. As 
a result, the complete thermal (and nonthermal) history 
of the universe in the context of the inflationary universe 
scenario may be dramatically different from the standard 
lore of the hot Big Bang cosmology. 
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APPENDIX 

Here we show how one can derive Eq. ([m|) or ( fj(i|) for 
the characteristic exponent fj,k from the analytic solution 
(pl|). For simplicity, we will consider here the case g 2 = A. 
Eq. ( |3l| ) describes both solutions, Xi(z) and X 2 (z). The 
resonant solution X(z) consists of four monotonic parts 
within a single period of the inflaton oscillation, see Fig. 
^. It turns out that at different quarters of the period 
either Xi(z) or X 2 (z) correspond to the exponentially 
growing solution. Indeed, the square of the resonant so- 
lution within the first quarter of a period is 



X 2 1/4 (z) 



Xl exp 



dz 



C 



M 1 (z)\ ^/z{l-z 2 ) 



(73) 



where Mi{z) is given by (29), C is given by (|32|), and X 2 
is the square of the resonant solution in the beginning of 
the period when z = 0. 

Within the second quarter one has 



X 2 /2 (z) = X 2 /4 (l) exp 



/J 



dz 



C 



Mi{z) \ ^/z(l-z 2 ) 



(74) 



Then the value of X 2 after half of a period is 

l 

Xl /2 {z = 0) = X 2 (z = 1) exp(-2C 



dz 



M 1 (z)^z(l~z 2 ) 
(75) 



where the integral is understood as its principal value. 
The resonant solution has the generic form X(z(x)) — 
P(z(x))e ,J - x , where P{z) is a periodic function. Since P 



18 



has a period equal to half of the period of the infiaton 
oscillation, Eq. (|75| ) is sufficient to find [i: 



= c f dz 

2 J M x {z)^z{l-z 2 ) 



> 0. 



The integral in this equation can be reduced to I(k 2 ) 
given by ©: 
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